Growth and Collapse of a Bose Condensate with Attractive Interactions 
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We consider the dynamics of a quantum degenerate trapped gas of 7 Li atoms. Because the atoms 
have a negative s-wave scattering length, a Bose condensate of 7 Li becomes mechanically unstable 
when the number of condensate atoms approaches a maximum value. We calculate the dynamics 
of the collapse that occurs when the unstable point is reached. In addition, we use the quantum 
Boltzmann equation to investigate the nonequilibrium kinetics of the atomic distribution during and 
after evaporative cooling. The condensate is found to undergo many cycles of growth and collapse 
before a stationary state is reached. 
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The recent observations of Bose-Einstein condensation 
in weakly interacting gases 0] has enabled a series of 
beautiful experiments that probe the dynamics of the 
condensate. The frequency and damping of the collec- 
tive modes of a condensate [0, propagation of sound in 
a condensate H , and recently, the growth of the conden- 
sate [[| , have been reported. Although these experiments 
have used atoms with positive s-wave scattering lengths, 
we show in this paper that the dynamical behavior of a 
negative scattering length gas, such as 7 Li, is especially 
interesting, and offers the opportunity to directly observe 
and study macroscopic quantum tunneling. 

A negative scattering length a implies effectively at- 
tractive interactions. In a spatially homogeneous gas, 
these interactions lead to ordinary classical condensation 
into a liquid or solid, preventing Bose condensation in 
the metastable region of the phase diagram However, 
confinement in an atom trap produces stabilizing forces 
that enable the formation of a metastable Bose conden- 
sate, if the number of condensed atoms is less than some 
maximum number N m . For a harmonic trapping poten- 
tial, Ruprecht et al. || showed that in mean-field theory 
N m ~ 0.57Z/|a|, where I = (Ti/mui) 1 / 2 is the extent of the 
one-particle ground state in the harmonic trap 0. For 
the 7 Li experiments of Ref. ||, N m ~ 1400 atoms, which 
agrees with the measured value. 

Although a condensate can exist in a trapped gas, it 
is predicted to be metastable and to decay by quantum 
or thermal fluctuations The condensate has only 

one unstable collective mode, which in the case of an 
isotropic trap corresponds to the breathing mode jl2],|l3| . 
The condensate therefore collapses as a whole, either by 
thermal excitation over a macroscopic energy barrier, or 
by quantum mechanical tunneling through it. The prob- 
ability of forming small, dense clusters is greatly sup- 
pressed because of the large energy barrier for this pro- 
cess, compared to that for the breathing mode. This 
suppression can also be understood from the fact that 
the typical length scale for fluctuations of the conden- 
sate is the healing length, which is precisely of the order 



of the condensate size near the instability point. 

Experimentally, it is also important to understand how 
such a condensate can be formed from a noncondensed 
thermal cloud by means of evaporative cooling. This 
question was recently addressed by Gardiner et al. in 
the context of the experiments with gases having a > 
fLU . These authors neglect the coherent dynamics of the 
condensate and focus instead on the kinetics of conden- 
sation (ll|. By treating the noncondensed atoms as a 
static 'heat bath' with a chemical potential that is larger 
than the ground state energy, they are able to derive a 
simple equation for the growth of the number of conden- 
sate particles that appears to fit well with experimental 
results [B. The same approach, however, does not work 
in the case of atomic 7 Li because it does not allow for 
the collapse of the condensate. Moreover, the process 
of evaporative cooling leads to dynamical changes in the 
noncondensed 'heat bath'. The study of both these ef- 
fects on the dynamics of the condensate is the main topic 
of this Letter. Some preliminary results have already ap- 
peared in a recent review article [[l6| . 

When collisions can be neglected, the dynamics of 
the condensate wavefunction ip for a gas with a < is 
well-described by the nonlinear Schrodinger (or Gross- 
Pitaevskii |L7|) equation 
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Here m denotes the mass of 7 Li, T 2B (0, 0; 0) = Anah 2 /m 
is the two-body T-matrix, and a is the scattering length 
of — 27 ao [fq|. For the external trapping potential V(x) 
we take a harmonic potential with an effective isotropic 
level splitting fiui — KfaxWyWz) 1 / 3 of 7 nK Note that 
we ignore the mean-field contribution from the noncon- 
densed atoms in this equation p9[ , because it is nearly 
constant over the size of the condensate and therefore 
only slightly affects the condensate dynamics. 

This description is semiclassical, and also neglects 
quantum and thermal fluctuations. These fluctuations 
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are most easily incorporated by means of the partition 
function of the condensate, which is a functional inte- 
gral / ef[^*]d[V>] exp(— S[ip*,ip]/K) over all periodic con- 
figurations of the condensate, with a weight determined 
by the (Euclidean) action of the nonlinear Schrodinger 
equation. It is most convenient to calculate this parti- 
tion function in terms of the density and phase of the 
condensate, defined by ip — J~pe %x . The Gaussian inte- 
gral over the phase field \ can then be performed exactly 
leaving only the determination of the functional integral 
Jd[p] exp(-%]/&) ; 

Unfortunately, this integral cannot be calculated ex- 
actly. However, since we are primarily interested in the 
dynamics of the unstable breathing mode of the conden- 
sate, we can use a variational method 1 11 2(J and consider 
only Gaussian density profiles 

p( X]q (t))=N (^ 7 -] 3/2 cx P (-^-] . (2) 
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Substituting this profile in the action S[p], we find that 
the dynamics of the condensate is equivalent to the dy- 
namics of a fictitious particle with mass m* — 3Nom/2 
in the unstable potential [|T) 



U(q) = N 



( zh 2 



\4mq 2 



3 2 2 

-rriLU q 
4 



iV h 2 \a 
\p2m mq 3 



(3) 



The rate of decay for both quantum tunneling and 
thermal fluctuations can now easily be calculated within 
this formalism Jll]] and are shown in Fig. [j]. For large 
numbers of condensate atoms, these collective decay 
mechanisms are much faster than the decay caused by 
inelastic two and three-body collisions, since the energy 
barrier out of the metastable minimum vanishes as Nq 
approaches N m . For experimentally relevant tempera- 
tures of 300 to 500 nK, the collective decay processes 
dominate for Nq greater than about 1100 atoms. 



Besides the calculation of the decay rates, the above 
anology also allows study of the dynamics of the collapse 
that occurs after the condensate has been driven out of 
its metastable minimum. A typical example is shown 
in Fig. ||. From this figure we see that the condensate 
first collapses with increasing speed along the potential 
hill outside the barrier. However, during the collapse, 
the condensate density grows rapidly, thereby increas- 
ing the decay rate for inelastic two and three-body colli- 
sions. Atoms that inelastically collide acquire subtantial 
energy and are ejected from the trap. Because of these 
loss mechanisms, the collapse is arrested when the num- 
ber of condensate atoms is of the order of one. Atoms are 
lost so quickly that the density of the condensate always 
obeys n|a| 3 <C 1, which is a necessary requirement for 
the validity of the nonlinear Schrodinger equation ([!]). If 
there were no inelastic collisions, the condensate would 
fully collapse J23|, and a more complex theory would be 
needed. 




FIG. 2. Typical evolution of the condensate during col- 
lapse. 
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FIG. 1. Decay rate of the condensate as a function of the 
number of condensate particles at 0, 100, 200, 300, 400, and 
500nK. The dashed line shows the decay due to inelastic two 
and three-body collisions |23] . 



The above remarks pertain to the dynamics of the con- 
densate in the absence of collisions with noncondensed 
atoms, but so far all experiments producing BEC have 
relied on evaporative cooling, which requires such colli- 
sions. Therefore, to investigate the kinetics of condensa- 
tion in a trapped gas, we use the Boltzmann transport 
equation in a way similar to the treatment of evapora- 
tive cooling in Refs. ^4j,^5|. In this approach, we de- 
fine a semiclassical distribution function /(x, p), with the 
number of atoms at position x with momentum p being 
dxdp /(x, p)/(27r?i) 3 . The evolution of the distribution 
function is then given by (2(| 

|/(x,p,t) = /(x,p)-r(x)/(x,p) , (4) 

where the derivative on the left is the total time deriva- 
tive. The effect of elastic collisions is given by 
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x <5 3 ( Pl + p 2 - p 3 - Pa) S(pI +P 2 2 -vl- Pi) 

x [/3/4(l + h){\ + h) - /i/ 2 (l + / 3 )(1 + /*)], (5) 

where the factors (1 + fi) account for the Bose statis- 
tics of the atoms. Here /, stands for /(x, pi) and 
a = 87ra 2 . Inelastic collisions, which lead to a loss of 
atoms from the trap, are described by T(x) = G\ + 
Gi J dp /(x, p)/(27rfi.) 3 , where the terms with G\ and 
Gi reflect the loss of atoms due, respectively, to collisions 
with hot background gas atoms and to inelastic collisions 
between trapped atoms. 

As for the case of a classical gas |^J,^5|, Eq. (^) can 
be simplified by taking the motion of the atoms to be 
ergodic, so that an atom with a given energy will sam- 
ple each available cell in {x, p} phase space with equal 
probability. The distribution function /(x, p) can then 
be expressed in terms of the energy distribution function 
f(E), through the relation 



/(x,p) = J dES(H(x,p)-E)f(E) 
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where H(x, p) = p 2 /2m + V(x) is the classical Hamil- 
tonian. The differential equation for f(E,t) is derived 
by substituting (||) into (||), as is described in detail in 
Refs. [ p4|j25| . The only difference here is the use of Bose 
statistics, which requires the insertion of factors l + f(Ei) 
as in Eq. (^). Although the semiclassical approximation 
results in a continuous function f(E), the quantum na- 
ture of the trapped gas can be recovered by restricting 
f(E) to values E n = (n + 3/2)tuv. 

The above model neglects the effect of the mean-field 
interaction energy. Per atom, however, this interaction 
energy is limited to about fous by the stability criterion 
arising from a being negative, which limits the effect of 
the mean-field energy on the kinetics of the gas. The 
accuracy of the approximations can be checked by con- 
sidering collisions between condensate atoms, since the 
mean-field interaction is largest in the condensate and 
the semiclassical approximation is least accurate for the 
ground state of the trap. We compare the rates for inelas- 
tic collisions between condensate atoms, G2 Jdxp(x) 2 , 
which depend only on the density. The exact condensate 
density is determined from Eq. (^) , using the variational 
method of Eq. (J2J) . In contrast, the semiclassical density 
derived from Eq. (j^) is 
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From these distributions, the ratio of the exact and 
approximate collision rates is 1.2(1 /q) 3 . When interac- 
tions are small, then q ~ I, and the error caused by the 
semiclassical approximation is only 20%. However, the 



mean-field interaction causes q to decrease as Nq is al- 
lowed to grow. From the variational calculation we find 
that q = 0.67Z just before the collapse. The error in the 
collision rate is then a factor of 4, and collision terms for 
other low-lying states will be incorrect by similar, though 
smaller, amounts. Although these errors are significant, 
they only occur for No relatively close to N m . For in- 
stance at Nq ~ 0.8 N m , the error is a factor of 2. As will 
be seen below, the condensate number is greater than 
0.8 N m for only a small fraction of time, so the approx- 
imations are useful, but for A^ — N m the quantitative 
predictions of the model will be somewhat inaccurate. 

The possibility for the condensate to collapse is in- 
cluded in the model using the decay rates given in Fig. [I]. 
A random number is chosen to determine whether a de- 
cay occurs during an integration time step, and when one 
does the condensate number is set to zero reflecting the 
rapid loss shown in Fig. |^. Also, evaporative cooling is 
included by setting f(E) = for E > E c (t), where E c (t) 
is the experimentally settable cutoff energy. 




FIG. 3. Typical evolution of condensate number No in re- 
sponse to evaporative cooling. The gas initially consists of 
4 x 10 5 atoms at 500 nK. During the first 5 seconds of evo- 
lution, roughly 40% of the atoms are removed by evaporative 
cooling, after which evaporation is halted. 

The response of the condensate number to evaporative 
cooling is shown in Fig. [| In this case, E c (t) consists of a 
5-second period of rapid cooling, after which evaporation 
is halted and the gas allowed to equilibrate. The subse- 
quent behavior shows the condensate alternately filling 
and collapsing, until finally the phase space density of 
the gas is too low to reach Nq — N m . The slow final 
decay is due to the trap losses. 

As shown in the figure, the time for which the oscil- 
lations in Ao persist is anomalously long compared to 
the elastic collision time r c = l/(nav) ~ 0.8 s. In or- 
der to investigate this phenomenon further, the density 
profiles generated by the model were fit to equilibrium 
Bose-Einstein distributions at various times. The val- 
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ues of x 2 resulting from the fits are plotted in Fig. ||. 
After reaching degeneracy, the curve shows an approxi- 
mately exponential approach to equilibrium, with a time 
constant of about 10 t c . In contrast, the same test per- 
formed on a nondegenerate cloud initially prepared in a 
nonequilibrium state yields an equilibration time of 5t c . 
The relatively slow approach to equilibrium for the de- 
generate case is presumably caused by the small phase- 
space volume of the condensate and the fact that due 
to the limit No < N m , the stimulated Bose scattering 
factors in Eq. (^) cannot become as enormous as they 
might in the a > case. Furthermore, the oscillations 
in No will persist until the distribution is very close to 
equilibrium, since the total number of atoms in the trap 
is much greater than N m . 
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FIG. 4. Equilibration of degenerate gas. Using the same 
conditions as in Fig. [], at each time the density distribution 
is fit to an equilibrium density distribution, and the resulting 
value of x 2 plotted. 

Note, finally, that the evolution in Fig. || represents 
only one possible outcome of evaporative cooling, and 
that because of the stochastic nature of the collapse a 
given evolution is not repeatable. Experimental measure- 
ments of No in the degenerate regime should, however, 
exhibit large fluctuations between and N m . Observa- 
tion of such fluctuations and measurement of their statis- 
tics would provide confirmation of the collective nature 
of the collapse. If sufficiently low temperatures can be 
reached, these fluctuations would be evidence of macro- 
scopic quantum tunneling. We therefore believe that a 
dilute Bose gas with a < presents unique opportunities 
for studying the dynamical properties of a condensate. 
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